A self-organized model for cell-differentiation based on variations of molecular decay rates

Systemic properties of living cells are the result of molecular dynamics governed by so-called genetic regulatory networks (GRN). These networks capture all possible features of cells and are responsible for the immense levels of adaptation characteristic to living systems. At any point in time only small subsets of these networks are active. Any active subset of the GRN leads to the expression of particular sets of molecules (expression modes). The subsets of active networks change over time, leading to the observed complex dynamics of expression patterns. Understanding of this dynamics becomes increasingly important in systems biology and medicine. While the importance of transcription rates and catalytic interactions has been widely recognized in modeling genetic regulatory systems, the understanding of the role of degradation of biochemical agents (mRNA, protein) in regulatory dynamics remains limited. Recent experimental data suggests that there exists a functional relation between mRNA and protein decay rates and expression modes. In this paper we propose a model for the dynamics of successions of sequences of active subnetworks of the GRN. The model is able to reproduce key characteristics of molecular dynamics, including homeostasis, multi-stability, periodic dynamics, alternating activity, differentiability, and self-organized critical dynamics. Moreover the model allows to naturally understand the mechanism behind the relation between decay rates and expression modes. The model explains recent experimental observations that decay-rates (or turnovers) vary between differentiated tissue-classes at a general systemic level and highlights the role of intracellular decay rate control mechanisms in cell differentiation.Comment: 16 pages, 5 figure

Complex and unexpected dynamics in simple genetic regulatory networks

Peer reviewedPublisher PD

The emergence of waves in random discrete systems

Essential criteria for the emergence of wave-like manifestations occurring in an entirely discrete system are identified using a simple model for the movement of particles through a network. The dynamics are entirely stochastic and memoryless involving a birth-death-migration process. The requirements are that the network should have at least three nodes, that migration should have a directional bias, and that the particle dynamics have a non-local dependence. Well defined bifurcations mark transitions between amorphous, wave-like and collapsed states with an intermittent regime between the latter two

Modeling Stochasticity and Variability in Gene Regulatory Networks

Modeling stochasticity in gene regulatory networks is an important and complex problem in molecular systems biology. To elucidate intrinsic noise, several modeling strategies such as the Gillespie algorithm have been used successfully. This paper contributes an approach as an alternative to these classical settings. Within the discrete paradigm, where genes, proteins, and other molecular components of gene regulatory networks are modeled as discrete variables and are assigned as logical rules describing their regulation through interactions with other components. Stochasticity is modeled at the biological function level under the assumption that even if the expression levels of the input nodes of an update rule guarantee activation or degradation there is a probability that the process will not occur due to stochastic effects. This approach allows a finer analysis of discrete models and provides a natural setup for cell population simulations to study cell-to-cell variability. We applied our methods to two of the most studied regulatory networks, the outcome of lambda phage infection of bacteria and the p53-mdm2 complex.Comment: 23 pages, 8 figure

Emergent Oscillations in Networks of Stochastic Spiking Neurons

Networks of neurons produce diverse patterns of oscillations, arising from the network's global properties, the propensity of individual neurons to oscillate, or a mixture of the two. Here we describe noisy limit cycles and quasi-cycles, two related mechanisms underlying emergent oscillations in neuronal networks whose individual components, stochastic spiking neurons, do not themselves oscillate. Both mechanisms are shown to produce gamma band oscillations at the population level while individual neurons fire at a rate much lower than the population frequency. Spike trains in a network undergoing noisy limit cycles display a preferred period which is not found in the case of quasi-cycles, due to the even faster decay of phase information in quasi-cycles. These oscillations persist in sparsely connected networks, and variation of the network's connectivity results in variation of the oscillation frequency. A network of such neurons behaves as a stochastic perturbation of the deterministic Wilson-Cowan equations, and the network undergoes noisy limit cycles or quasi-cycles depending on whether these have limit cycles or a weakly stable focus. These mechanisms provide a new perspective on the emergence of rhythmic firing in neural networks, showing the coexistence of population-level oscillations with very irregular individual spike trains in a simple and general framework

Probability distributed time delays: integrating spatial effects into temporal models

Background: In order to provide insights into the complex biochemical processes inside a cell, modelling approaches must find a balance between achieving an adequate representation of the physical phenomena and keeping the associated computational cost within reasonable limits. This issue is particularly stressed when spatial inhomogeneities have a significant effect on system's behaviour. In such cases, a spatially-resolved stochastic method can better portray the biological reality, but the corresponding computer simulations can in turn be prohibitively expensive.Results: We present a method that incorporates spatial information by means of tailored, probability distributed time-delays. These distributions can be directly obtained by single in silico or a suitable set of in vitro experiments and are subsequently fed into a delay stochastic simulation algorithm (DSSA), achieving a good compromise between computational costs and a much more accurate representation of spatial processes such as molecular diffusion and translocation between cell compartments. Additionally, we present a novel alternative approach based on delay differential equations (DDE) that can be used in scenarios of high molecular concentrations and low noise propagation.Conclusions: Our proposed methodologies accurately capture and incorporate certain spatial processes into temporal stochastic and deterministic simulations, increasing their accuracy at low computational costs. This is of particular importance given that time spans of cellular processes are generally larger (possibly by several orders of magnitude) than those achievable by current spatially-resolved stochastic simulators. Hence, our methodology allows users to explore cellular scenarios under the effects of diffusion and stochasticity in time spans that were, until now, simply unfeasible. Our methodologies are supported by theoretical considerations on the different modelling regimes, i.e. spatial vs. delay-temporal, as indicated by the corresponding Master Equations and presented elsewhere

Avalanches in a Stochastic Model of Spiking Neurons

Neuronal avalanches are a form of spontaneous activity widely observed in cortical slices and other types of nervous tissue, both in vivo and in vitro. They are characterized by irregular, isolated population bursts when many neurons fire together, where the number of spikes per burst obeys a power law distribution. We simulate, using the Gillespie algorithm, a model of neuronal avalanches based on stochastic single neurons. The network consists of excitatory and inhibitory neurons, first with all-to-all connectivity and later with random sparse connectivity. Analyzing our model using the system size expansion, we show that the model obeys the standard Wilson-Cowan equations for large network sizes ( neurons). When excitation and inhibition are closely balanced, networks of thousands of neurons exhibit irregular synchronous activity, including the characteristic power law distribution of avalanche size. We show that these avalanches are due to the balanced network having weakly stable functionally feedforward dynamics, which amplifies some small fluctuations into the large population bursts. Balanced networks are thought to underlie a variety of observed network behaviours and have useful computational properties, such as responding quickly to changes in input. Thus, the appearance of avalanches in such functionally feedforward networks indicates that avalanches may be a simple consequence of a widely present network structure, when neuron dynamics are noisy. An important implication is that a network need not be “critical” for the production of avalanches, so experimentally observed power laws in burst size may be a signature of noisy functionally feedforward structure rather than of, for example, self-organized criticality

MicroRNA-Mediated Positive Feedback Loop and Optimized Bistable Switch in a Cancer Network Involving miR-17-92

MicroRNAs (miRNAs) are small, noncoding RNAs that play an important role in many key biological processes, including development, cell differentiation, the cell cycle and apoptosis, as central post-transcriptional regulators of gene expression. Recent studies have shown that miRNAs can act as oncogenes and tumor suppressors depending on the context. The present work focuses on the physiological significance of miRNAs and their role in regulating the switching behavior. We illustrate an abstract model of the Myc/E2F/miR-17-92 network presented by Aguda et al. (2008), which is composed of coupling between the E2F/Myc positive feedback loops and the E2F/Myc/miR-17-92 negative feedback loop. By systematically analyzing the network in close association with plausible experimental parameters, we show that, in the presence of miRNAs, the system bistability emerges from the system, with a bistable switch and a one-way switch presented by Aguda et al. instead of a single one-way switch. Moreover, the miRNAs can optimize the switching process. The model produces a diverse array of response-signal behaviors in response to various potential regulating scenarios. The model predicts that this transition exists, one from cell death or the cancerous phenotype directly to cell quiescence, due to the existence of miRNAs. It was also found that the network involving miR-17-92 exhibits high noise sensitivity due to a positive feedback loop and also maintains resistance to noise from a negative feedback loop

Delay-Induced Transient Increase and Heterogeneity in Gene Expression in Negatively Auto-Regulated Gene Circuits

A generic feature in all intracellular biochemical processes is the time required to complete the whole sequence of reactions to yield any observable quantity-from gene expression to circadian rhythms. This widespread phenomenon points towards the importance of time delay in biological functions. Theoretically time delay is known to be the source of instability, and has been attributed to lead to oscillations or transient dynamics in several biological functions. Negative feedback loops, common in biochemical pathways, have been shown to provide stability and withstand considerable variations and random perturbations of biochemical parameters. The interaction of these two opposing factors-of instability and homeostasis-are features that are widespread in intracellular processes. To test the effect of these divergent forces in the dynamics of gene expression, we have designed and constructed simple negatively auto-regulated gene circuits consisting of a basic regulator and transcriptional repressor module, and compared it with one, which has delayed repression. We show, both theoretically and experimentally, that delayed repression induces transient increase and heterogeneity in gene expression before the gain of stability effected by the negative feedback. This design, therefore, seems to be suitable for conferring both stability and variability in cells required for adaptive response to a noisy environment